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I. INTRODUCTION 

Recently Georgi raised the interesting possibility that physics beyond the Standard Model (SM) may contain 

a sector that is conformally invariant in the IR region (guaranteed by a zero of the beta function), and classically 
scale-invariant in the UV; we refer to these as the unparticle iU) and Banks-Zaks {BZ) phases, respectively. The 
transition region between the two phases is characterized by the scale of dimensional transmutation Ku A specific 
realization of this idea can be found in [J] ; following this reference we will assume that the new sector is described as 
an asymptotically free gauge theory in the BZ phase. 

This novel idea has received substantial attention within the high-energy community, mainly in connection with 
the phenomenology of such models. Here we discuss some fundamental issues in the evolution of the Universe in the 
presence of this type of new physics (though studies of the cosmological consequences of the proposal have appeared in 
the literature 0-0, these publications ignore several essential aspects which are discussed below). In sec.|ll]we derive 
an approximate equation of state for the NP sector. Then, in sec. Illll we use this together with the expected SM-NP 
interactions P, to determine the conditions under which the SM and NP sectors were in equilibrium. In sec. IIVI 
using the experimental constraints derived from Big-Bang Nucleosynthesis (BBN) we obtain non-trivial bounds on 
the parameters of the theory. The Appendices A and B are devoted to presentation of two alternative derivation of 
the Boltzmann equation. 



II. THERMODYNAMICS OF UNPARTICLES 



In order to understand the thermodynamic behavior of the new sector ^ we use the expression for the trace anomaly 
of the energy momentum tensor of a gauge theory where all the renormalized masses vanish [TT| : 



0^ = Yg^[FrFa,A , (1) 



'Electronic address: bohdan. grzadkowski @fuw.edu. pi] 
^Electronic address: jose.wudka@ucr.edu 
^ A similar idea was discussed also in Q . 

The thermodynamics of conformal theories has been studied extensively [To| . but these results have been apparently ignored where 

unparticles are concerned. 
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where (3 denotes the beta function for the couphng g and A'^ stands for the normal product. 

The basic assumption for the unparticle phase is that the (3 function has a non-trivial IR fixed point at 5 = 7^ 0. 
Modeling the unparticle sector by a gauge theory, we assume that for low temperatures ^ 

/3 = a((?-(?0, a>0, (2) 

in which case the running coupling reads 

gitJi) = + un"; f3[g{p)] = au^i" , (3) 

where u is an integration constant and ^ is the renormalization scale. 

We look for the lowest-order corrections to the conformal limit (where 6*^ = 0) when the system is in thermal 
equilibrium at temperature T, is isotropic and homogeneous, and does not have any net conserved charge. Since (3 
vanishes in the conformal limit, in ^ we can take (TV [F^'''Fa ^i/]) equal to its conformal value (we denote the thermal 
average by (• • ■)); taking the renormalization scale ^ ~ T wc then expect 

{N [FrFa ^.]) = bT'+'> , (4) 

where 7 is the anomalous dimension of the operator. Using /^^^ = Pu ~ 3P;^, where pu and Pu denote the energy 
density and pressure of the unparticle phase, together with ([3)) and (j4]) then gives 

Pu-3Pu^AT^+'- (^A^^^, S^a + ^y (5) 

where we took /i = T. 

Combining ([5]) with the thermodynamic relation d{pV) + P dV = T d(sV) (s is the entropy density), when p and 



P are functions of T only'', and integrating, we find 



4+S 



Pu = laT^+ i]T^+' 



su = \aT^ + A{l + -\T^+' (6) 



3 \5 

where a is an integration constant and we assumed 5 ^ Q. 

It is worth noticing that the terms cx A correspond to deviations from the standard relativistic relation V oc T~^. 
The behavior at low temperatures depends on the sign of 5, we will assume 5 > 0. Then 



"iPu = Pu 



1 - ^Pu 



exhibiting the lowest-order corrections to the often-used expression P ~ wp, w =const. This effect might be of interest 
in the discussion of the possible dark-energy effects contained in this model, but will not be discussed here. 

Elucidating the cosmological effects of the modified equation of state ^ lies beyond the scope of the present paper, 
we merely remark that the NP increases the coefficient of the T'^ term in p and induces 0{T^) corrections; e.g. in the 
radiation-dominated era the scale parameter behaves as (1 -I- cT^Y^'^ /T (c =const.). 

In general we expect A oc since Ay is the scale associated with broken scale invariance; then the energy density 
for the new sector in the unparticle phase equals 

Pu = —T 



T « Au (8) 



where we replaced a = S^ir/tt^ (hereafter we use the normalization from Maxwell-Boltzmann statistics) and giR , the 
effective number of relativistic degrees of freedom (RDF) , will be estimated below. 



^ The cases where /3 has a higher-order zero at can be treated similarly. 
* A consequence of having assumed the absence of net charges. 
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In the BZ phase we assume the theory is asymptotically free so that, up to logarithmic corrections, 

PBz = ^9BzT^; T » Au (9) 

where gsz denotes the RDF in this phase. 

For intermediate temperatures the explicit form of the thermodynamic functions requires a complete non- 
perturbative calculation and the choice of a specific model; fortunately we will not need to consider the detailed 
behavior of the system. Given that p oc in both the IR and UV regions, for our purposes it will be sufficient to 
use the interpolation 

PNP = ^,9Npr^ gNP = gBz0{T - Au) + gue{Au - T) (10) 

where gu — [giR + (T/Au)^ f] while NP stands for 'new physics'; ^np will be continuous at T = Au when / = 
gsz — giR, which we now assume. It is worth noting that a mass distribution of unparticlcs with the spectral density 
cx; [ij generates the term oc / in (fTO|) with S = 2{du — 1), assuming that the contributions with fi > T 

decouple. We emphasize that pUj) will be used only as a rough but convenient approximation that reproduces the 
expected behavior at low and high temperatures. In cases of interest we expect giR ^ g^z 3* f so that the terms 
oc are subdominant. 

Estimating gi^ directly form the model Lagrangian is a non-trivial exercise, due to the expected strong-coupling 
nature of the theory in the infrared. Using, however, the AdS-CFT correspondence [l^l we find 

5IR = ^{LMpif (11) 

where L denotes the AdS radius of curvature and Mpi is the Planck mass. Given that L is expected [13] to be 
significantly smaller than 1/Mpi, it is justified to expect that 

.giR>O(100) (12) 

In the following we will use this as our estimate for the RDF in the unparticle phase. 

In order to estimateog^ one must specify the details of the non-Abelian theory in the ultraviolet regime. For the 
models considered in [J] we find 

9BZ - 100 (13) 

This result is based on a model for which the couping constant stays within the perturbative regime throughout its 
evolution. There is also non-perturbative lattice evidence [l3| that gauge theories exhibiting an infrared fixed point 
obey (fT3|l . In the following we will adopt this estimate. 

The energy density pn was also discussed in , however the expression presented in this reference agrees with ([8]) 
only when giR = and therefore does not include the leading low-temperature behavior of the theory. 



III. SM-NP INTERACTIONS AND EQUILIBRIUM 



The presence of a NP sector of the type considered here can have important cosmological consequences since, even 
when weakly coupled to the SM, its energy density will affect the expansion rate of the universe; this can then be used 
to obtain useful limits on the effective number of degrees of freedom ^np- This calculation requires a determination 
of the relationship between the temperature of the NP and SM sectors to which we now turn. 

The interactions we will consider have the generic form 

Ant = cOsmOnp (14) 

where the first term is a gauge invariant operator composed of SM fields (possible Lorentz indices have been sup- 
pressed), while the second operator is either composed of BZ fields or is an unparticle operator, depending on the 
relevant phase of the NP sector. The coupling e in general has dimensions and is assumed to be small. For the specific 
calculations presented below we will assume for simplicity that Osm.np are both scalar operators. 

Leading interactions involve SM operators that can generate 2 particle states since states with higher particle 
number will be phase-space suppressed. From such interactions we obtain the NP<->SM reaction rate P, which will be 
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precisely defined below. The two sectors will then be in equilibrium whenever r>_ff, where H denotes the Hubble 
parameter [l^ . and decouple at the transition temperature Ty: 



r2 Stt 



T = Tr. T^H; H'^—^p,,,- (15) 



where 



Ptot = PSM + Pnp 
3 

= ^ (ssmTIm + 5npTnp) (16) 

We denote by Tsm and Tnp the temperatures for the SM and NP sectors which can be different when these sectors 
are not in equilibrium 

The approach to equilibrium can be described using either the Kubo formalism (appendix [X| or a suitable extension 
of the Boltzmann equation formalism (appendix [B| . It follows form the expressions derived in the appendices that 
the conditions near equilibrium are determined by the equation 



i) + AHi)^~T&, 1?- Tnp -Tsm (17) 



where, using the Kubo formalism, 



^ " iS^ + ''^|/ "^^i dt^d^x ^OgM(-is,x)e)sM(t,0))(oNp(-«s,x)ONp(t,0))| (18) 

The Boltzmann equation (BE) calculation also yields p?)) with the rate given by 

^ = {— + — ) ^ E / d^^^d<^SMP{EsM - E'^^fe-f"^-^ \M\^ {2^f5{KsM - K^p) (19) 

12T-' V5SM 5NP/ 2T ^_^J 

where Ai is the matrix element (with no spin averaging) derived form the SM-NP interaction Lagrangian ^, E^m 
and E'o^yi denote the initial and final energies of the Standard Model particles in the reaction, and -STsm.np the total 
4-momenta of each sector for the reaction; we have also assumed the Boltzmann approximation (neglecting Pauli 
blocking or Bose-Einstein enhancement) and denoted by d^'sM.NP the appropriate phase-space measures (without any 
spin factors). In particular, for the unparticle phase we use [l| 

d^U = A,„e(q°) 9{q^) (g2)<i»-2_A_ (20) 

where A„ = (47r)3-2"/[2r(n)r(n - 1)]. We show in appendix |B] that (HHD and HH) are, in fact, equal. 

The solutions to pT|) yields p ex R~'^ in the absence of the collision term (proportional to F), as expected for a scale 
invariant theory. It is also important to note that, in contrast to other authors (0-0)1 CHI) contains an unparticle- 
decay term (see appendix [B)) . as we find the arguments (based on the deconstruction picture [10) for neglecting these 
contributions unjustified^. 

The detailed calculation of F requires a specific form of the interaction OsmOnp (see above for a specific example). 
However for the purposes of the remaining calculations only the basic properties of F, such as its dependence on T 
and the relevant RDF will be needed. These properties can be obtained using dimensional analysis: if the dimensions 
of the operators are, respectively dgM and c?np and if the number of degrees of freedom involved in this interaction 
are g'^y^ and g'^p, then, including a phase-space factor we find 

2 \ f I 

p ^ ^ ^fftot y2dsM+2dNP-7. ^ gSM ^NP i^^X 

^ (47r)"SM+nNp-l ' - gg]^ gjjp ' 



^ The SM-SM and NP-NP interactions are not included because of our assumption that each sector is in equiUbrium: these processes are 
much faster than the ones generated by I I14I I and insure that each sector has a well-defined temperature at all times. 

® II19I I gives the same result within the unparticle scenario or the deconstruction approach; in the latter case the vanishingly small 
coupling constant of the deconstructed field is compensated by the large number of particles of the same invariant mass in the initial 
state. Unparticle decay was discussed recently in [l6l . 
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where nsM and tt-np denote numbers of SM and NP fields in the corresponding operators; in the unparticle phase we 
take g'^p = and tt-np = '^{du — l)j where du denotes the dimension of Ojj. 

The value of A depends on the details of the model. Above the Higgs (0) mass (we assume ~ w = (0)) the 
most important operator is Osm = (f>^(f>', in this case g'^^ = 4, so A ~ (4/(7sm) • (.9np/.9np)- Below there are many 
dimension 4 SM operators relevant for the SM-NP equilibration, e.g. £4)e (containing an extra suppression by the 
factor ^ v/Mu; (,e denote a lepton isodoublet and isosinglet respectively), or B^^B^" (where B is the hypercharge 
gauge field), in this case we expect g'gj^ ~ gsM, so that A ~ 5np/5np- 



A. The Banks-Zaks phase. 



We will assume that the BZ sector corresponds to an SU (uc) Yang-Mills theory with Uf vector-like massless fermions 
in the fundamental representation (denoted by qsz). Assuming that Ku > v, the leading SM^NP interaction is of 
the form 

= ^ ('/'V) {qBzQBz) (22) 
where we assume that all flavors in the BZ sector couple with the same strength. In this case (e = 1/Mu) 

Denoting by T^z f the solution to (fT5|) when F is given by and imposing also the consistency conditions 

Mu > Tisz-f > ^u, we obtain ( gtot is evaluated at T^z-f) 



Mu Agtot Mpi Mu 



B. The unparticle phase. 
In this case we will consider only interactions of the form [H (fc = dsM + dsz — 4) 



^ = ^^77^0smOu (25) 



Using (PT|) we obtain (here we use tinp = 2{du — 1)) 



\ A / \ \ 2k / rri \ 2dsM+2du~7 



(47r)nsM+2d„-3 y ^yY, 



(26) 



Denoting by Ty^f the solution to pSj) when P is given by ((26)) . and imposing also the consistency condition Ku > Tu_f, 
we obtain (here gtot is evaluated at Ti(_f) 



Ti 



u-f 



(47r)"sM-H2d„-3 Mu^^^ 



< 1 (27) 



For du < 4.5 — dsM, T/ H has the singular property of increasing as T drops, whence SM and NP will equilibrate 
for T < Tu-f (thaw- in); due to the constraints ^ on djj {du < 1 is excluded [o]) this can only happen for Osm = 4'^(l>- 
The opposite occurs if du > 4.5 — dsM (freeze-out). For du = 4.5 — dsM, the approximations (fT6|l . ([2T|l are insufficient 
and a detailed calculation is required to determine freeze-out and/or thaw- in conditions; we will not consider this 
special case further. 

There are various possible scenarios for decoupling of the NP sector. The situation in the very early Universe 
(T > Mu) depends on the UV completion (including the mediator interactions) of the NP and will not be considered 



^ The bounds on strictly hold in the conformal limit; we expect deviations oc g{T) — g» ~ (T/A./j)"' which we neglect. 
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FIG. 1: Regions in the Au — Mu plane corresponding to various freeze-out and thaw-in scenarios for du ~ 3/2, 2, 3, 7/2. 
Dark grey: SM-NP decouphng in the unparticle phase only; light grey: no SM-NP decoupling; in the white regions Tu_f < v 
[Au, Mu are in TeV units). We assumed gsyi = gBz = gu = 100, psM = 4, g'sz ~ 50 and glj — du. For the BZ phase: 
nsAi = nNP = 2, dsM ~ 2 and djvp = 3, while for the U phase: nsAi = 2, njvp — 2{du — 1), dgm = 2 and djvp = du- 



here. If (|24p holds then we have a standard freeze-out scenario: the SM and NP sectors will be in equilibrium down 
to T ^ Tisz-f and decouple below this value; thereafter the two sectors evolve keeping their entropies separately 
conserved. Since no mass thresholds or phase transitions are crossed * the SM and NP temperatures remain equal 
down to T ^ Au ■ 

The situation for Ky^T is more complicated. If (P7)) holds (which defines a region in the Ay — Mu plane), 
decoupling occurs in the unparticle phase. For T > v the most relevant operator is Osm — <i'^<t>: and both thaw-in 
(for du < 2.5) and freeze-out (for du > 2.5) may be present. For w > T all the relevant SM operators have dsu — 4, 
and only freeze-out is possible; in this case Tu-f may be significantly smaller than v. 

Other parameter values lead to more complicated scenarios, e.g. a double decoupling: freeze-out in the BZ, thaw-in 
in the unparticle phase and then freeze out below v. In spite of the many possibilities, there is always a temperature 
below which the SM and NP decouple. 

In Fig. [1] we show regions in the (Au,Mu) space that correspond to various freeze-out and thaw-in scenarios for 
a reasonable parameter choice. For this calculation we assumed that Osm ~ <i>^ <t' is responsible for maintaining 
the equilibrium between the SM and NP (so dsM = 2). For consistency that choice implied an additional constraint 
Tu-f > V (below V other SM operators arc relevant). For interactions with the BZ phase an operator oc {(f)^ 4)){qBZ(lBz), 
was adopted (in which case d^z = 3). 



We neglect the possibility of right-handed neutrino decoupling. 
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IV. BIG BANG NUCLEOSYNTHESIS 



The light-element abundances resulting from BBN are sensitive to the expansion rate that determines the temper- 
ature of the universe (see e.g. [11]), which can be used to restrict possible additional RDF, or, in our case, gi^. We 
express our results in terms of the number of extra neutrino species, ATV^, defined through 



3 7/4 



^^P^^il^nj'^'^^'^^'' ^^^^ 

annihilation (T^ stands for the photon temperature). For AN^ we adopt the 
recent bounds obtained in [18|: AA^^ ~ 0.0 ± 0.3atat(2cr) ± O.Bgyst- 

We first consider the case where SM and NP were in equilibrium down to a temperature Tf > v, and decoupled 
thereafter. Then the entropy conservation for the NP and SM sectors implies 

9^p{Tf){TfRf)' = g*MTNp){T^pR)' 

gsuiTfmRff = g*sM{T^){T^Rf (29) 

where Rf is the scale factor at the decoupling while R corresponds to temperature of photons (T^p is the cor- 
responding NP temperature); g^p and g^-^^ stand for the NP and SM effective numbers of RDF conventionally 
adopted for the entropy density. After e~^e^ annihilation neutrinos and photons generate the dominant SM contribu- 
tion, but their temperatures differ. Using standard expressions we find 

9sm{T^) = 37 \ ' 30) 

g-r + ge 

where gi stands for the number of RDF corresponding to the species i. Assuming that gNp is almost constant in 
the temperature range we are interested in and neglecting possible right-handed neutrino decoupling effects, the two 
sectors had the same temperature down to the electroweak phase transition; thereafter the temperatures split as the 
SM crossed its various mass thresholds and the entropy was pumped into remaining species. Entropy conservation 
([29)1 in both sectors then implies 



Pnp — 



g-i gii.e.vy^'^ 



gj + ge gsM{v) 



(31) 



where gsuil, e, = g^ + ge + gi^, while gsuiv) stands for the total number of SM RDF active above T ~ v. Note that 
the above relation holds regardless if the decoupling happened during the BZ or unparticle phase. Then combining 
with (1^51) we obtain 
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gsM{v) 



gsM{l,e,v) 



AN, (32) 



Using the standard expressions for the SM quantities [Tj| the BBN constraint on AN^, then implies gm ^ 20 at 95% 
CL. It is worth mentioning here that F measures the decay rate of unparticles into SM states. After decoupling, when 
T < H these decays become very rare (the NP — > SM life-time becomes larger than the age of the universe ^ ^^/H). 

More severe constraints could be obtained if NP and SM remained in equilibrium down to the BBN temperature. 
That occurs for Ky, My ~ TeV and du ~ 1; the relevant operator being Bf^^B^^Ou- Then, since temperatures of 
the NP and SM sectors are the same, one obtains 

'?i« = lM^ ^^'^ (33) 

4 Kg-r + geJ 

which leads to 5ir<0.25 at 95% CL. 

When decoupling occurs between v and Tbbn the bound on gm lies between 0.25 and 20. When the SM and 
NP are never in equilibrium the BBN constraints can be used to bound pnp, but not gm since Tnp is then not 
known. These bounds should be compared to giR> 100 typical of specific models [J e.g. for an SU{3) gauge theory 
with 16 fundamental fermion multiplets, and expected from AdS/CFT correspondence [13]. We conclude that many 
unparticle models will have difficulties accounting for the observed light-element abundances. 
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V. SUMMARY 



Using the trace anomaly we argue for a form of the equation of state for unparticlcs that contains power-hke 
corrections to the expression for relativistic matter; this allows us to determine temperature dependence of the 
energy and entropy density for unparticles. We then derive the Boltzmann equation for the BZ phase and postulate 
a plausible form for this equation for unparticles; using this we determine the conditions for NP-SM equilibrium. 
Finally we derive useful constrains on the NP effective number of degrees of freedom imposed by the BBN. 
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APPENDIX A: DERIVATION OF THE REACTION RATE USING THE KUBO FORMALISM 



In this section we follow closely the arguments presented in [l^ . We consider a thermodynamic system, not necessar- 
ily in equilibrium, with macroscopic observables {ai} associated with operators {ci}. We assume the thermodynamics 
of the system is described by a density matrix p 



p = exp 



where the pi and /3 are parameters, and = ri(/i, /3) is a function chosen such that tvp = 1, that is 
The pLi are determined by the condition 

ai = Tvpa^ = - -— 
\opt 



(Al) 



(A2) 



(A3) 



It is important to note that p differs from the usual grand-canonical density operator in that the are not assumed 
to be conserved, so the ai will not be constant: 



a^{t) = Tr{p a^{t)} = Tr{p(i) a,}; a^{t) = e'"*a,e-'"' 



iHt „JHt 



, p{t)=e-^"'pe 



ai{t) denotes the average of at time t for a distribution for which the average of at t = is ai = ai(0). 
We now assume the pi are small, then a straightforward calculation yields 



(A4) 



(A5) 



where, for any operator ^, 



Now let 



a'^{t) = Trpa^(t) = Tr {pe'"'a'^e-'"'} ; a', = a, ~ (a,) 



(A6) 
(A7) 
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so that, to first order in /i, 

4it)^J2 f ds{a'^{-is)ar{t))p., (A8) 



Using now the cychc property of the trace, {(^{z)ri{z')) = (^(z — z') rj) — ri{z' — z)) for any operators 77 and any 
complex times z, z' . From this it follows that 

^ {a[{-^s)a'^{t)) ^ - {a[{~^s^{t)) , (A9) 



hence 



(AlO) 



Next, using the definition 

m^t[H,m (All) 

and the cyclic property of the trace, 

ds {a[i~is)d'j) = -I ([a-, a;-]) = -i {[a,,aj]) . (A12) 



Collecting all results and using a'^ ~ hi, 



ij IJ-j 



J 

= / d-"^ I dt.il \ {aj{-is)ai{t)) + i {[at,aj]) 



which is the celebrated Kubo equation. It is important to note that the r ^ limit is subtle }19j . 
Suppose that the system is composed of two sub-systems, labeled '1' and '2' with a Hamiltonian 



(A13) 



H = Hi + H2 + eH' ■ [Hi,H-2]^Q, e<l (A14) 

and take ai = Hi , 02 = H2', in this case p describes two systems at different temperatures that weakly interact 
through eW . Then 

a, - {H,) = Vp, (A15) 

where denotes the energy density and V the space volume of the system. We imagine that each subsystem has a 
well defined temperature Ti but that these change slowly due to the presence of H'\ we also require the systems to 
be close to equilibrium with each other so that |T — T^J <C T. In this case the left hand side of (|A13p corresponds to 
a'^ while on the right hand side we can take the t — > cxd limit since the integrand is damped at times larger than the 
characteristic times of systems 1 and 2; see Ref. [l^ for details. In this case 

= c,5T, ; 5T,^T,-T (A16) 

where Ci denote the heat capacities per unit volume at temperature T . 
When e = 0, the density matrix (|A1|) becomes 

p\ _p = g/30-/3(l-,ii)ifi-/3(l-M2)ff2 (^A17) 

which corresponds to non- interacting subsystems at temperatures T!j = T/(l — /i^), whence 

Z^^^^-JT,, T=i; (e = 0) (A18) 
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Then ((IT3l) gives 

Vc,6T, = ^-Y_^g,j6T,- g,,^J dsj dt(^H,{-is)H,{t)) (A19) 

where 

H, = t[H,H,]^te[H',H,] ^ = e'^^ij.e-'^^ = 0(e) , (A20) 

so that Q is of order e^; since we work to the lowest non-trivial order in H' , this also justifies the use of (jAlSp . 
Now we need to evaluate Q. Following (|14p. we assume 



H' = - j d^xCiCa (A21) 

then 

i[H',Hi],=^ = j d^^i[Hi,0i]02^ j d'i^di02 (A22) 
and, similarly, i[H,H2] — J d^x020i ^From this 

{ ^ (HiHs)Hi{t))^ = J d^x (e)i(-j,s, x)di{t, y)) {02{~is, ^)02{t, y)) 
|^(i/i(-zs)ij2W)} = / d^^d^y ldi{-is,^)Oi{t,y)) l02{-is,^)d2{t,y)] 
I ^ (ff2(-*s)iji W)} - y ^^'x d-V x)Oi(t, y)) ld2{-is, x)02(t, y)) 

{^(^2(-»s)ij2(t))| = I d^xd^y {0,{~is,^)0,{t,y)) (^d2{-ts,x)d2{t,y)) (A23) 

where the (• • ■) separates into a product because when e = averages separate into averages over systems 1 and 2 
which are independent. For the case where the Oi are scalars and even under time reversal all the above correlators 
are equal up to a sign, so that 

g = ^^Gv(\ (A24) 



where V denotes the volume of space and 

G::^J^dsJ^ dt J d^x (^Oi{-is,x)di{t,0)'j (^O2{-is,x)d2{t,0)'j (A25) 

Substituting ((J24l) in (MM gives ciSTi = -028X2 = -{e^G/T){5Ti - ST2), then 

dt{5T,-ST2) = -r{5T,-ST2); T =(- + -) (A26) 

\ci C2J T 

The quantity G can be evaluated using the tools of finite-temperature field theory. To facilitate this let 

Jo = -iOi dt O2 (A27) 
then, setting e = and using invariance under space translations. 



1 



/3 







G ^ -^J ds I dt I d3x(Jo(-is,x)Jo(t,0)) 



= l^^j"^ ds J d^xdte~''^'^''^''^e{t){.M-is,0)Joit,x))\ (A28) 
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In this form G can be evaluated in terms of the correlator of two Jq currents. We took the real part, which is the one 
that yields the relevant width, and introduced fc as a regulating 4-momentum. The limit w, k — > requires care, for 
the present case one should first set k = and then take lo to zero [llj . 

In order to compare this result to the one derived using the Boltzmann equation it proves convenient to do a 
Lehmann expansion of G, which involves matrix elements of the form Jo|m). In terms of Feynman graphs, such 
matrix elements will include pieces that are not connected to Jq; these disconnected pieces factorize and cancel the 
factor exp(/351o) 2l| that appears in the definition of the average (|A6p . We then find 



/5|(^ko|m)_|^(27r)4<5W(p„-p„0 



(A29) 



Up to now we have assumed that the volume of the system is kept fixed, but this can be easily relaxed. The 
calculation involves obtaining the thermodynamic potential to order /i^ and will not be presented here, the final result 
is the expected one: the time evolution equation becomes pi + AHpi = — ^Qij5Tj/T where V /V ~ ~3H. 



APPENDIX B: THE BOLTZMANN EQUATION 



We again imagine two sectors, labeled 1 and 2; within each the interactions are strong enough to maintain equilib- 
rium at temperatures Ti 2; the sectors interact only though p4p . We denote by the distributions of particles a 
in sector i; the corresponding Boltzmann equation is 



' [dx- J ^^^^^ 



CM 



(Bl) 



where the right hand side denotes the collision term. 

We consider first a process of the form Xi + X2 X[+ X2, where Xi, X[ [i = 1, 2) denote states in system i. If 
a particle labeled by a^^^ is in Xi, then the corresponding collision term C[i^q(i)] is given by 

K(i)] = - j d<^'x,x\ \M{X ^ X' f {2^)H{K^ +K2- K[ - K!,Wx,x' 



Cx,X'['^a(-i-'> 

Mx.x 



K' 



n (li'^&d)) n (i±^c(2)) n ^d^^^ n ^<^^'^ 

bWex[ c(2)ex^ ci(i)eXi e(2)eX2 

n (li'^bd)) n (i±^c(2)) n n ^^cd 

bWeXi c(2)eX2 dWex[ eWex!^ 

(ii;,,K,), E,^ ^a<'' ^ =Y. ^-^'^ 

{E^K'^, Ei= J2 ^ai^> ^ k,(., 



(B2) 



where d^'x x' denotes the corresponding invariant phase space measure for all particles except a^^^ (as indicated by 
the prime), the Lorentz-invariant matrix element, and Ea, k^ denote the energy an momentum of particle a. The 
upper sign corresponds to bosons, the lower to fermions. 

We will assume spatial homogeneity, so that the v will depend only on time and energy, and also assume kinetic 
equilibrium, so that the density functions take the usual Fermi-Dirac or Bosc-Einstcin form, but with time dependent 
temperature and, possibly, chemical potential. Then 



Mx,x' 
vx,x' 



^-Ei/Ti-E2/T2 _ g 



-E'JTi-E'2/T2 



vx.X'] 

n (i±^6(^)) n (i±^c<i)); 



(B3) 



Using this we can derive the time dependence of the energy density; for simplicity we will carry out the calculation 
in flat space. The energy density associated with the a(i) is 



(2^ 



EaWVaW - 2 / d<^a{l)E^^^)V^^l) 



(B4) 
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Integrating (jBl|) over p we find 

dtPai^)\x,x' = d^x,x'E^,^■, \M{X X'f {2^)H{K^+K2-K[-K'^)Afx,x' , rf^x.x' - d^'^^^'d^ai^) ; (B5) 

where the notation on the left hand side indicates that this corresponds to the change in p^(i) generated by this 
particular X ^ X' reaction. The total time derivative is obtained by summing over all states X, X' such that 
a(i) G Xi: 

Pan-> = - J2 I d-i>x.X'E,a) \MiX ^ X')f {2^)H{Ki + K2 - K[ - K^)Mx^x' (B6) 

X,X'; (a(i)eXi) 

The time derivative of the total energy density for each sector is then obtained by now summing over all a^-'^' : 

Pi = -J2 f d'^x^x'Ei \M{X ^ X')f {2^f5{Ki + - K[ - K^)AAx,a-' (B7) 



To make this look more symmetric consider the contribution with X and X' exchanged. Since \M.\^ is the same 
but TV changes sign we can write 

Pi = -\Y1 f d^x.x'{Ei-E[)\M{X^X'f{2Tif5{Ki+K2-K[-K'2)Nx.x' (B8) 

X,X' 

The corresponding expression for p2 is obtained by switching the 1 and 2 subscripts. 

We are interested in cases where the Maxwell-Boltzmann statistics are adequate, so vx,x' — 1, and when the 
temperatures are similar: Ti = T + 6Ti. Using the energy conservation condition Ei + E2 = E[ + E2, we find 

Mx^x' ^ -e-(^^+^=)/^^i^ {6T, - 5T2) (B9) 
Also, ignoring non-relativistic contributions to the energy density 

P^ = cM (BIO) 

where q is the heat capacity per unit volume. Collecting all expressions gives 

dt {5Ti -5x2) = -r (dTi -ST2) , 



r = 



- -1 - 

ci + C2 / 2T 



f d^x,x'P{E^ - E[fe-^^- \M{X ^ X'f {2^)H{Kx - Kx') ■ (Bll) 

X',X 

In order to compare this with the Kubo formula we use 

MiX^X') = {X'\Cint\X)^e{X'\(Di02\X) (B12) 
where we work to lowest non-trivial order in the interaction. Using Jo , defined in (|A27p , we find 

(^Vol^).=o - 2(£;i - E[) {X'\0,02\X)^^, (B13) 
where we took e = since we are interested only in the leading contributions to T. Then 

,Cl C2, 



X',X 

Using then the Lehmann expansion ()A29p we find 



J2 I d^x,x'e~^''- K^Vol^)!' i27rfSiKx - Kx') (B14) 

Vt V ^ 



Cl C2 T 
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exactly as in the Kubo formalism . 

Despite its intuitive appeal the Boltzmann approach contains conceptual difhculties for the case of strongly inter- 
acting theories, for which concepts such as the particle densities Va are ill defined. In this case the definition of T 
(|A26[) obtained through the Kubo equation is preferable where the relevant matrix elements can, at least in principle, 
be obtained numerically. 
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We have used the Boltzmann approximation in identifying En in I IA29I I. which is the total energy of state |n), with Ex which is the 
sum of the energies of the particles in state \X). These energies are approximately equal for a sparse system, where this approximation 
holds. 



